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A first order ■perturbation expansion is carried out in order to analyze the 
effect of small spatially uniform strains on the lowest order {even) TE and 
TM modes in an anisotropic dielectric waveguide. This generalizes the 
results of an earlier paper in which the effects of certain special cases of 
uniform strain were calculated. Unlike in these special cases, the perturbed 
modes are, in general, neither purely TE or TM, and one effect of two of the 
offdiagonal components of the strain is to tilt the plane of polarization and 
change the relative phase of the two polarizations. To first order, the modes 
are not exponentially attenuated. Some numerical examples are considered 
in order to illustrate the results. It is found that, under appropriate condi- 
tions, the effect of the small strain may be quite large in relation to its 
magnitude. 

I. INTRODUCTION - 

The concept of a multilayered dielectric waveguide is central to 
the theory of the GaP electro-optic diode modulator. 1 - 9 As part of a 
detailed study of the properties of electro-optic diode modulators, 
Nelson and McKenna 4 have investigated the possible discrete modes 
which can propagate in a number of such waveguides and have cal- 
culated the detailed properties of the lowest order mode of each 
polarization. 

In the fabrication of a p-n junction a certain amount of strain is 
always introduced. Because of the photoelastic effect 10 this strain will 
induce a change in the dielectric matrix describing the unstrained p-n 
junction. In general the strain will be spatially nonuniform, making 
it extremely difficult to calculate modes in such a structure. However, 
a knowledge of the effect of a spatially uniform strain on the mode 
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structure would provide insight into the effects of nonuniform strain. 
The effects of certain special cases of uniform strain on the modes of 
a simple model of a dielectric waveguide were calculated in Ref. 4. 
In the present paper we complete this investigation and calculate the 
modes in the same model dielectric waveguide when subjected to an 
arbitrary uniform strain. We use first order perturbation theory in a 
small parameter describing the magnitude of the strain. 

Although the work presented in this paper was motivated by re- 
search on the theory of the electro-optic diode modulator, the results 
have considerable relevance to the theoiy of the GaAs injection laser. 
Here too, various dielectric waveguide models have been used to 
explain the light containment. 11-14 The same problems of strain exist, 
and the results of this paper give a qualitative picture of the effects 
of strain on modal structure. The effect of strain on completely 
different types of electro-optic light modulators have been studied by 
Kaminow 15 and by DiDomenico and Anderson. 16 

II. FORMULATION OF THE PROBLEM, AND RESULTS 

In Ref. 4 the symmetric step model was used to study the effects 
of strain. This very simple model exhibits many of the main features 
of interest in dielectric waveguide models. 

The model consists of an anisotropic crystalline slab bounded by the 
planes x = ± w, whose refractive index is raised uniformly by some 
constant amount, the physical origin of which is still obscure, and which 
is embedded in an isotropic medium of relatively lower index of refrac- 
tion. The central slab represents the junction region whose anisotropy is 
caused by the junction field Ej acting through the electro-optic effect. 9 
The direction of the rc-axis is always taken parallel to E/. The isotropic 
medium represents the normal GaP. 

The model is determined by its dielectric matrix, which in the absence 
of strain and for certain orientations of E y with respect to the crystal 
axes can be diagonal in a coordinate system having its z-axis parallel to 
"Ej. For such orientations of Ej, the diagonal matrix elements of the 
dielectric matrix K^ix), a = 1, 2, 3, depend only on x. (We use x, y, z 
for the coordinates rather than x,, x 2 , x 3 .) The matrix elements in the 
absence of strain are then defined by the equations 

K[°\x) = K a , | x | < w (1) 

K™(x) =K , | x | > w (2) 

where a = 1, 2, 3 (see Fig. 1). 



WAVEGUIDE STRAIN 



1935 



K ( „ 0) (x) 





(a) 



(b) 



Fig. 1 — (a) The coordinate system used in the symmetric step model, (b) A 
graph of K a (x). 

There are two orientations of E y of particular interest which allow us 
to diagonalize the dielectric matrix in the desired coordinate system. 
If Ej is in the [111] direction, then the x, y, and z axes can be taken in 
the [111], [IlO], and [Tl2] directions, while if Ejis in the [100] direction, 
the x, y, and z axes can be taken in the [100], [Oil], and [Oil] directions. 
The set of axes determined by the unstrained model will be used in all 
the strain calculations and the dielectric matrix will always be referred 
to these axes. See Ref. 4 for further details of the model. 

In the presence of a uniform strain, the dielectric matrix is in 
general no longer diagonal, and we can write for the dielectric matrix 
elements K a p (x) , 



K aa (x) = KT(x) + v S aa 



a = 1,2, 3 



(3) 



K a ,(x) = V S aP , **($. (4) 

The symmetric matrix {t]S a p) is the contribution of the photoelastic 
effect 10 which we have written in this form for convenience in the 
perturbation analysis. The matrix elements S a p are spatially constant. 
We assume that n~ 2 S a p is of order unity, where n is the index of refrac- 
tion of GaP and ?/ is a small parameter. In Section II we express rjS a p 
in terms of the strain matrix and give estimates for the size of -q. 
We now seek solutions of the Maxwell curl equations 



V X E - -m„H, 



(5) 
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V XH = e K(x)-E, (6) 

of the form 

E = e(x) exp z'(o/ - pkz), (7) 

H = h(.r) exp i(pt - 0kz). (8) 

These solutions correspond to waves travelling in the positive z direc- 
tion, where k = wO^o)' = 27r/X is the free space wave number and X the 
free space wavelength of the light. 

In the strain free cases (77 = 0) , there are both TE and TM modes 
and these modes can be either even or odd functions of x. At most 
only a finite number of modes can exist, and Ref. 4 shows that for the 
typical parameter values encountered in GaP diode modulators only 
the lowest order even TE and TM modes can exist. For that reason 
we confine ourselves here to solutions which in the limit of zero 
strain (77 = 0) reduce to even modes. However, the perturbation 
technique used here applies equally well to solutions which in the 
limit 77 = reduce to odd modes. 

When 7} ?± 0, we seek solutions of Maxwell's equations of the form 

e a (x) = A a exp —kp(x — w) 

+ B a exp —kq(x — w), x ^ w (9) 

= C a exp kr(x + w) + D a exp ks(x -\- w), x ^ —w (10) 
= F a exp ikfx + G a exp —ikgx 

+ L a exp iklx + M a exp —ikmx, \ x \ ^ w (11) 

for a = 1, 2, 3. The general solution in each region is a sum of four 
linearly independent solutions, but in the regions \x\ > w, the bound- 
ary conditions at infinity eliminate two of these solutions. The expres- 
sions for h a (x) can be obtained from equation (5). The various coef- 
ficients and parameters A a . . . , p, . . . can be expanded in powers of -q 

A a = AT + tjX" + ••• , (12) 

P = Vo + vPx + • • • . and so on. (13) 

In Section III we list the terms in these expansions of order zero and 
one in 77, and in Section IV we outline their derivation. In this section 
we merely discuss some of the features of the solutions. 

We refer to solutions which in the limit as 77 — > reduce to even 
TE modes as "perturbed TE modes"; similarly, we refer to "perturbed 
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TM modes." Expressions for the unperturbed TE and TM modes 
(rj = 0) are given in equations (25) through (32). 

Although the expressions we have obtained for the coefficients and 
parameters are quite complicated, several features of the perturbed 
modes stand out. If <S 12 ^ or £ 23 5^ 0, then the modes cannot be purely 
TE or TM. In general an important effect of the strain is to tilt the plane 
of polarization. This tilt is in general a function of x but not of z. Since 
the coefficients A a , ... are complex in general, the relative phase 
between E y in the perturbed TE mode and E z in the perturbed TM mode 
is a function of x. This relative phase at x = cannot be determined 
unless the method of excitation is known, since all the A a cannot be 
determined, as is shown in Section IV. Ref. 4 considers the special case 
where <S 12 = S 2Z = 0, S 13 ^ and shows that the modes are rigorously 
TE or TM. That paper calculates only such parameters as /3 and p, not 
such coefficients as 4 a and B a . The parameters are expanded in two 
small quantities 8 and A describing the unstrained dielectric matrix. If 
we expand the expressions for the parameters in this paper to first order 
in the same small quantities 5 and A (to second order for 0) complete 
agreement is obtained with the Ref. 4 results. 

In the absence of strain, the surfaces of constant phase for both 
TE and TM modes are the planes z = constant. However, in the 
presence of strain, the surfaces of constant phase are no longer planes, 
and are different for the perturbed TE and TM modes. 4 

Finally, j3 + ypi is real in all cases. Thus at least to first order in 
?/ the modes experience no exponential attenuation as they propagate. 
In order to get some feel for the magnitude of the effects involved, 
we consider several numerical examples. We first estimate the order 
of magnitude of 77 by relating it to observable phase differences. 
Consider a plane wave whose free space wavelength is A propagating 
over a distance I in a medium of index of refraction n + An. The phase 
difference Ap which this wave would experience over the same wave if 
the index of refraction were n is 

A<p = ^Z(An). (14) 

A 

If 77ft. 2 is the photoelastic contribution to the dielectric constant, then 

A?i = (n a + ijw 2 ) 1 — n = \tit]. (15) 

Therefore, we have 

r? = \A<p/nlir. (16) 
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Now the upper limit of phase shifts observed 17 in GaP at X = 6328 A 
over a length I = 0.6 mm is about t/4. Taking 18 n = 3.31 this yields 
7] ~ 0.8 X 10~\ This probably represents an extreme upper limit, and 
so we will assume ij = 10~ 6 in our examples. Recent X-ray measure- 
ments 10 of strain in P doped Si yield a value of -q of about 10~ 6 when the 
concentration of the dope is N D tt 10 18 atoms per cubic centimeter. 9 

It can be shown that the matrix (7jS a g) is approximately related to the 
strain matrix (e aj3 ) by the equations 10 



vS aP = -n* £ P, 



(3/n.Cj, 



(17) 



where n is the index of refraction and P aPliV are the elasto-optical coeffi- 
cients. Crystals of class 43m have only three different elasto-optic 
coefficients when referred to the crystal axes. (See p. 251 of Ref. 10.) 
For GaP these are 20 

P„ = -0.151, P 44 = -0.074, P 12 = -0.082. (18) 

Since 18 n = 3.31 for GaP, and n~-S a p is at most of order one, it 
follows that the magnitude of the strain is roughly proportional to 77. 
In order to obtain the values of the elasto-optic coefficients in the 
coordinate system used in this paper, it is necessary to make a trans- 
formation of the elasto-optic tensor from its representation in the 
crystal axes. We will not do that here; rather we take (rjS a p) as given. 
In Table I we define three possible strain contributions to the dielectric 
matrix, labelled a, b and c. Matrices a and b were chosen to demon- 
strate the effect of the off-diagonal elements Si 2 and <S 2 3, respectively 
(Ref. 4 considered the effect of S±s alone) , while c was chosen to dem- 
onstrate a possible effect when all the off-diagonal elements are nonzero. 

For a GaP diode modulator we can write 2, 4 

K a = n 2 (l - 8 a ), a = 0, 1, 2, 3, (19) 

where n = 3.31 is the index of refraction of GaP 17 and the quantities 
8 a , a = 1, 2, 3 are functions of the applied bias voltage V. In the original 



Table I — Strain Contribution to the Dielectric Matrix* 



Type 


V 


Sn 


Sa 


S33 


Sn 


Sa 


•Sis 


a 


10- 8 











10 








b 


10" 8 














10 





c 


lO" 6 


7.07 


7.07 


7.07 


7.07 


7.07 


7.07 



Components of the strain contribution, <S,y , and the magnitude parameter t). 
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symmetric step model 8 is independent of V. For Ej \\ [111] Ref. 4 
showed that 

5, = -25, 5 2 = <5 3 = 5, (20) 

where 

5 = - Eyr 41 n 2 /(3)% (21) 

.and where E/ is the (spatial) average junction field, r tl is the electro-optic 
coefficient, and n is the index of refraction. For a typical diode (diode 
KC46CA of Ref. 9), Ej (measured in V/m) is related to the diode half 
width w (measured in m) and bias voltage 7 (measured in V) by 9 

E, = (2 " V)/(2w). (22) 

The half width can be determined by capacitance measurements 8 - n 
and related to the bias voltage by 

w(V) - 0.139 X HT 6 (1 - 7/1.8) - 380 . (23) 

Using the value /- 41 = -0.86 X 10" 1 - m/V,* we can now calculate 
8i , 8 2 , and 8 :t as functions of V. 

For this diode, 8 = 1.612 x 10 :1 . However, it has been shown 
that the voltage dependence of the parameters of the symmetric step 
model is not correct, and the double walled waveguide much more 
closely describes the true voltage dependence. 4 ' ° We have not used 
the double walled model because it is analytically complex. Instead, 
since the modes in the single and double walled guides are very similar 
in form because they both decay exponentially as functions of x out- 
side the guide, we have used the single walled model but simulated 
the voltage dependence of the double walled model. This has been 
achieved by letting 8 vary with voltage. The voltage variation of 8 
has been obtained by requiring the equality of expressions (2.33) and 
(3.18) in Ref. 4 for the decay constants p, and letting W\ = itf(0) and 
io 2 = iv {V) . This yields the relation 

Su - (2.24 X \QT w )/w. (24) 

In Table II we list these basic constants describing the unstrained 
diode as functions of V. Using these values, we can calculate from 
equations (33) through (37) the parameters of the unstrained TE 



* This is the undamped value of r« given in Ref. 18. After these calculations 
were made it was determined that the clamped value r« = —0.97 X 10" 12 m/V 
should be used. However, since our results supply only qualitative information 
about actual diodes, we have not redone the numerical example. 
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Table II — Characteristics for a Typical GaP Diode* 



Bias 

voltage 
(10 


to 
(10-* cm) 


10< «i 


10* 6: 


10* it 


10« 5. 


-2 
-12 
-24 


1.87 
3.04 
3.85 


-1.17 
-2.51 
-3.67 


0.58 

1.25 

1.84 


0.58 
1.25 
1.84 


12.00 
7.38 
5.82 



* Given as functions of the applied reverse bias voltage V. The half width of the 
junction is w, and the components of the unstrained dielectric matrix are Kj = 
n 2 (l — Sj),j = 0, 1, 2, 3, where n = 3.31 is the index of refraction of GaP. 

and TM modes for a = 6328 A. These values are listed in Table III. 
Finally, in Tables IV and V we list the parameters of the correspond- 
ing perturbed TE and TM modes respectively. The accuracy of those 
terms less than 10 -4 is uncertain in case c of Tables IV and V. In 
Figs. 2 through 7 we plot some of the components of the perturbed 
TE modes correct to first order in rj. In Figs. 2, 3, 5, 6, and 7 the 
imaginaiy part of the component is negligible and is neglected, while 
in Fig. 4 the real part is negligible with respect to the imaginary part 
and is neglected. In all cases the e 3 component is negligible compared 
to the ei component. We have chosen the undetermined coefficients 
so that at x = 0, z = 0, e 2 in the perturbed TE mode and e x in the 
perturbed TM mode have zero phase to first order in 77. 

This example illustrates how much tilting of the plane of polariza- 
tion, or coupling of the TE and TM modes, is to be expected. The 
£12 component produces the main effect, which from Figs. 2 and 3, 
is a maximum tilt of the plane of polarization of 3.5°. This effect 



Table III — Unperturbed Mode Parameters* 



Type 


Bias 










voltage 










mode 


(V) 


/So 


PO 


/o 


h 


TE 


-2 


3.308 


0.0226 


0.1096 


0.1180t 


TE 


-12 


3.309 


0.0195 


0.0796 


0.1022+ 


TE 


-24 


3.309 


0.0160 


0.0641 


0.1007+ 


TM 


-2 


3.308 


0.0259 


0.1089f 


0.1174 


TM 


-12 


3.309 


0.0307 


0.0759f 


0.0994 


TM 


-24 


3.309 


0.0363 


0.0552t 


0.0953 



* Describing the unstrained TE modes, and the parameters /3 , p and l de- 
scribing the unstrained TM modes as functions of the applied reverse bias voltage V. 
The wavelength of the light is 6328 A. 

t Derived parameters Jp for the TE modes and/o for the TM modes. These derived 
parameters appear only in first and higher order corrections to the field. 
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Table IV — Parameters for Perturbed TE Modes* 



Type 
of 

strain 


Bias 

voltage 
00 


/?■ 


pi 


9i 


h 




/re(.pil 


/m(pi) 


Re(qi) 


Jm(gi) 


mi 


a 


-2 
-12 
-24 








221 
257 
312 




-221 
-257 
-312 














b 


-2 
-12 
-24 










1.51 

1.51 

1.51 




-1.51 
-1.51 
-1.51 












c 


_2 
-12 

-24 


1.07 
1.07 
1.07 


156 

1S1 
221 


-0.3X10- 7 
+0.2X10-' 
-0.2X10-' 


-156 
-181 

-221 


-2.14 
-2.14 
-2.14 


2.13 
2.13 

2.13 


-0.4X10-' 

-0.4X10- 8 
-0.4X10- 8 



* For all perturbed TE modes fi = gi = 0. 

decreases with increasing reverse bias voltage. However, it should be 
noticed from Figs. 4 and 5 that the coupling effect resulting from 
&> 3 increases with reverse bias voltage. The e 1 component is roughly 
proportional to rj, so a doubling of the strain would double the mode 
coupling. Mathematically, the existence of this relatively large effect 
results from the largeness of the factor c given in equation (64) for 
perturbed TE modes and in equation (84) for perturbed TM modes. 
The TM modes exhibit a similar behavior. 

From Tables IV and V we see that the changes in the parameters, 
VVi> VQu 7 //?i> and so on, are indeed small, which gives us confidence 
that the perturbation treatment is reasonable. 



III. FORMULAS FOR THE SOLUTIONS 

To list the formulas for the coefficients and parameters, A a , . . . , 
p, . . . (which appear in the expressions (7) through (11) for the 
solutions in terms of the various parameters describing the symmetric 
step model and the strain matrix) , we begin by writing down the solu- 
tion for the strain free (?; = 0) case for both the even TE and TM 
modes. When ?j = 0, we have for the even TE modes 



ei(x) = e 3 (x) = 0, all x 

e 2 {x) = cos (Jcf Q x), I x I ^ w 

= cos (fc/ow) exp kp (w — \ x \), 
while for the even TM modes 



(25) 

(26) 

x I ^ w (27) 
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Fig. 2 — The relative amplitude of the real part of ei for the perturbed TE mode. 
e x {x) = cos (kl x) , | x | ^ w (28) 

a; |), | x | > w (29) 

(30) 
x I £ w (31) 



= -=r cos (/cZ„w) exp /kp («' 
e a (a;) =0, all x 



e :i (x) = i-^-rr s'm (kl x), 

= i V £f cos (kl w) sgii (a:) exp kp n (w - \x\), \x\^w. (32) 

Pfrfi-0 

The parameters in these equations are given for the TE modes by 
the positive roots of the system of equations for p , po, and / 

pi = &l — K , (33) 

f = K 2 - £ , (34) 

/„ tan(/ar/ ) = p (35) 

while for the TM modes the parameters are the positive roots of the 
system of equations for p , fio, and / , consisting of equation (33) and 

11 = K 3 - {KJKM (36) 
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Fig. 3 — The relative amplitude of the real part of ea for the perturbed TE mode. 
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Fig. 4 — Tho relative amplitude of the imaginary part of e t for the perturbed 
TE mode. 
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Fig. 5 — The relative amplitude of the real part of e s for the perturbed TE mode. 

K lo tan(AwZo) = K 3 p . (37) 

The expressions for h(x) can be obtained from equation (5). 

We now turn to listing the formulas for the coefficients and param- 
eters of the solutions for the perturbed TE and TM modes. For both 
the perturbed TE and TM modes we have the relations 

p» _ (&V - (-VV„(,s 22 + W + &(&* - W + 2*p /Ui 3 ] 

9i \po/ \4p„Ao/ 

-f iKoWoiSw + foA.] 1 } 1 , (38) 

where pi corresponds to the " + " sign and #! to the " — " sign. It is 
also true, at least to first order in rj, that 



e a (-x) = e a (x)*, a = 1, 2, 3, 



(39) 



hence we only list those parameters determining the solution in x ^ —to. 
For the perturbed TE modes Po, fo, and /? are the positive solutions 
of the system of equations (33) through (35), and Iq is then given 
in terms of /? as the positive root in equation (36). The remaining 
parameters are 



9o = Po , 
0o = /o , 



(40) 
(41) 
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Fig. 6 — The relative amplitude of the real part of e, for the perturbed TE mode. 



0, = S a2 /(20 a ). 



(42) 
(43) 



The quantities pi and q^ are now determined by equations (38) and 
(43) . Next we have 



/i - 0i - (£22 - 200/30/(2/0), 



(44) 
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Fig. 7 — The relative amplitude of the real part, of e 3 for the perturbed TE mode. 
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U = TTlh- [S U (K 3 - - 2K a fi p i + S, 3 (ff, - fi) ± 2&.M. 

(45) 

where I] corresponds to the " + " sign and m x to the "-" sign, and in 

equations (44) and (45) fii is given in equation (43). Notice that 
from (43) and (44) ,f 1 = g 1 = for the TE case. 
The expressions for the coefficients are 

A{ 0) = -/3 (/3 n £ 12 + ipoSn) cos(kUw)/[2K oVQ (pi - ffOl. ( 46 ) 

A< 0> = -q t cos (kfow)/( Pl - g,), (47) 

A 3 (0) = i(po/Po)A?\ (48) 

jj{°> = -A[°\ (49) 

B' ' = p, cos (kf w)/( Pl - <?,). ( 5 °) 

Bi 0) = i(po/Pn)B[°\ (51) 

Fl 0) = (?1 0) = 0, a- 1,3 (52) 

F vn = gm = i (53) 

Li 0) = Ml 0) =0, a - 1,2.3 (54) 

F{1) = Jc [=F SM - S V2 (K, - f )] , (55) 

or 

F 2 n = G?\ (56) 

Fi<1) = \c [=F S l2 o / O - &.(X, - /©] , (57) 
GT 

L *" = |J cos(Mo)laS 23 [2i3 2 (^ 1 - K )c - 1] 

± ^-° ^ I2 [2pS(^o - # 3 )c - 1]) , (58) 
Vo 

L\ l) = (PoK 3 /l K l )U 1 \ (59) 

Ml" = -{PoKJUKJM^ , (60) 

L' 1 ' = M' 1 ' = 0, (61) 
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where 

a = [l K cos (kwl ) + p K 3 sin (kwk)]~ l $ (62) 

b = [loKos'm (kwl n ) — p Q K 3 cos (kwl )]~ l , (63) 

c = (#,#3 -Kafi-K&y 1 , (64) 

and pi and q 1 are the values appropriate to the TE modes. 
Finally, we write down the three combinations 

A[ l) +B{" = -i{^/p )[A^ 4-53 1 '] 

- cos (kf w)[(pl - K )S 12 + ip p S 23 ]/(2K pl), (65) 

AP + B l 2 " = cos (ki w)[F^ + G 2 X) ], (66) 

A l a l) + Bj" = [F™ + GH cos (kf w) + iffiS" - Gj ,} ] sin (*/ u>) 

+ [L' 1 ' + M? } ] cos (M„w) + t[l4" - JfJ 1} ] sin (fcZ u>). (67) 

The coefficients Fi" = G"', and hence A"' + B ( 2 l) , are arbitrary and 
correspond to an overall multiplicative constant. They can be set equal 
to zero with no loss in generality. We discuss this point further in 
Section III. Moreover, the individual coefficients A, (1) , A^, B^, and 
B% 1) cannot be determined at this stage. However, the terms we have are 
sufficient to determine each component of the field up through order 
one in t\. 

For the perturbed TM modes p Q , h, and /?o are the positive solu- 
tions of the system of equations (33) , (36) , and (37) , and / is given 
in terms of /? as the positive root in equation (34). The parameters 
g and p are still related by equation (40), m and lo by equation 
(42), and g and / by equation (41). The remaining parameters are 

ft - [l/(2/3 )][tfo(#iZo 2 + K a pl)l+ f]" 1 

• {[Snti/KiWoK&l + K\ll + f] 

+ [K&X/KtW^Kt - K )pl + rll, (68) 
where 

f = (kp w)(K 2 3 pl + KUl). (69) 

With the aid of (68) and (69) for fi lt p x and q x are determined by 
equation (38) , /i and g x by equation (44) , and l x and m x by equation 
(45). 
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The expressions for the coefficients are 

A™ = -ft cos (kwl )(2p Pi ~ 2p p. - £ 22 )/[2poA'o(p, - ffl )], (70) 

A 2 (0) = -ft cos (kwl o )(0 o S VJ + ip o S 23 )/[2p o o Ko( Pi - ?,)], (71) 

B[ m = ft cos(kwl )(2p (3r - 2p g, - flW/PWfo(Pi - fc)L (72) 

J5<o> = _ A (o) ( (73) 

where {J 1 is given by equations (68) and (69), pi and <?i are the values 
appropriate to the TM modes, and A^ 0) is related to A[ 0) by equation 
(48) and £< 0) to B[ 0) by equation (51). Furthermore, 

F™ = Gl 0) = 0, a = 1,2,3 (74) 

LJ ' = Mr = I, (75) 

L<°> = Mi 0) = 0, (76) 

H m = -ilf< 0) = (Z ft)/(2/3 ft), (77) 

F{" = Fj" = (?r = G' 1 ' = 0, (78) 

Fi = cos(M ){,S a& 2 [ft + 2cpl(K - ft)] 

± pobSnlKt + 2c/3a 2 (A', - K )]}/(4p (3 K ), (79) 



L20> = -4c[5„ ± Z o ft£ 23 /03 o ft)] (80) 



Ai° + #" = [F, (n + G 2 (1) ] cos (fcw/ ) + t[F 2 u - <# n ] sin (fcw/o) 
+ [L 2 (1) + Mi 1 '] cos (fcw/ () ) + i[U l) - M^) sin (fctoZo), (81) 
where 

a = \p cos (kwj ) — /o sin (fcw/o)] -1 , (82) 

b = [/„ cos (kwj ) + p sin (kwfo)]' 1 , (83) 

c = (ft - ll - ffiT 1 . (84) 

Just as in the perturbed TE case, the coefficients cannot all be deter- 
mined uniquely. We can with no loss of generality set 

Z4 U = -AT*" = 0. (85) 
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Once this choice is made, we have 

L«> = [S 13 - ZoA - kp + (S 33 - 2l l l )(l K l /p K 3 )]/(2(3 l ), (86) 
M\ X) = [-5 13 - l^ - m,j8 + (£33 - 2l m 1 )(l K l /p K 3 )]/(2p l ), 

(87) 

+ iWkwiK^lSn/K, + KXS 33 /K 3 - 2K a p <t p l )]/(K a p )} cos (M ) 

(88) 
4<» + b»> = -tGWPb)[Aj 1} + #"] + (#i C08(kwlo)/K o PoPo) 
'[-KMpo - (SoSnCpo - K )/(2p K ) 
+ PoPoS 33 /(2K ) - iplSn/Ko]. (89) 

A knowledge of these terms is sufficient to determine each component 
of the field up through order one in 77. 

IV. DETAILS OF THE CALCULATIONS 

In standard fashion H can be eliminated from equations (5) and (6) 
by taking the curl of equation (5) and by making use of the assumed form 
of the solutions, equations (7) and (8). There results the system of 
equations 

ifi%*-f? ei + t,K la e a = 0, (90) 

# - $V. + t, K 2a e a = 0, (91) 

;S L + tf|f+ t^ a e a =0, (92) 

where we have introduced the new independent variable 

t = kx. (93) 

The standard boundary conditions 21 on E and H yield the conditions 
that e 2 , e z , (tei/dl-, and de 3 /d£ + ipei must be continuous at £ = 0- = kw. 
The general plan of the calculation is first to consider the equations 
obtained by substituting into equations (90) through (92) the expres- 
sions for e a in the various regions given by equations (9) through 
(11). From these equations, one can determine up through first order 
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in 77 all but one of the parameters and some of the coefficients as 
functions of the parameter ft. Upon substituting these values into 
the boundary condition equations, a set of equations is obtained from 
which ft and some of the remaining coefficients can be determined. 
Since the components of the electromagnetic field satisfy a linear, 
homogeneous system of equations, it follows that if e a (x) , a = 1, 2, 3 
is a solution set, then so is (1 + a x -q + a&f + . . .)e a {x), a = 1, 2, 3, 
where the constants a lf a 2 , . . . are arbitrary. For example, if the 
coefficients A a , B a , . . . given by expansions of the form (12) rep- 
resents a solution, then the coefficients given by expansions of the form 

A a = AT + vi^AT + AT] + • • • , (94) 

with the same a-i used in each expansion, represent another solution. 
Thus unless the corresponding zeroth order coefficient is zero, the first 
order coefficient cannot be uniquely determined. We do, however, 
have the arbitrary constant a x at our disposal. The multiplicative 
constant (1 + (hy + • • •) can onr y ^ e determined from a knowledge 
of the excitation of the mode. 

If the assumed expressions for e a in £ ^ a given by equation (9) are 
substituted into equations (90) through (92) we get the set of homoge- 
neous, linear equations in A a , a = 1, 2, 3, 

(K Q + r,Su - /3 2 M, + r,S 12 A 2 + (r,S ia - i(3p)A 3 = 0, (95) 

vSuAi + (K + r,^ + V 2 ~ /*V 2 + i& u Ai = 0, (96) 

( V S» - ifip)A x + r,S 23 A 2 + (K + v S-, 3 + p 2 )A 3 = 0, (97) 

plus a similar set of equations with A a replaced by B a and p replaced 
by q. The condition that these equations have a non trivial solution, 
the vanishing of the determinant of coefficients, yields a relation be- 
tween ft and p of the form 

D(p, 0) = 0, (98) 

where D (p, (3) is a quartic polynomial in p and ft. The second set of 
equations involving the B a and q yields the same determinantal equa- 
tion with p replaced by q, 

D(q, 13) = 0. (99) 

That is, q is a second root of the quartic. If p, q, and ft are expanded 
in powers of ?; as in equation (13), equations (98) and (99) can be 
expanded in powers of 7/ and the coefficients of the various powers of 
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7/ can be equated to zero. The vanishing of the lowest order term 
yields equation (33) , which is satisfied by both p and q , thus also 
yielding equation (40). The vanishing of the first order coefficients 
shows that p x and q^ are the two roots of a quadratic which are 
given by equation (38). These results are independent of the TE or 
TM character of the mode. 

Equations (95) through (97) can now be expanded in powers of 77 
by substituting in the expansions of p, ft, and A a . The three similar 
equations involving q, ft, and B a can be expanded in powers of 77 in 
the same way. Because p and ft satisfy equation (33) , equation (96) 
vanishes to zeroth order in 77, while equations (95) and (97) yield 

(K - pl)A< 0) - iftoVoAT = 0, (100) 

-i/3oPoA\ 0) + (K + pl)Ai 0) = 0. (101) 

The determinant of this pair of homogeneous equations vanishes because 
equation (33) is satisfied, so a nontrivial solution exists. The quantities 
B[ 0) and B ( 3 0) satisfy the same equations. Using equation (33), it follows 
from equation (101) that A\ 0) and A { 3 0) are related by equation (48), and 
Bl 0) and B 3 0) are related by equation (51). 

To first order in 77, equations (95) through (97) are 

(K - &Al li - iftoPoAl" 

= ~[(S U - 2ft (3 1 )Ai 0) 4- S l2 A^ 0) 4- (£ 13 - ipofti - ip,ft )A ( 3 0) ], 

(102) 

S 12 Al 0) + (2p oVl - 2^/3, + S 22 )A ( 2 0) + S 23 A ( 3 0) = 0, (103) 
-tfoMi" + (K 4- pl)Al" 

= -[(&, - tpoft - ipiPo)Al 0) + S 23 A< 0) 4- (S 33 + 2poP.)Aj 0) ]. 

(104) 

With the replacement of A a by B a and p by q in equations (102) 
through (104) we obtain the first order equations satisfied by the B a - 
At this stage we must differentiate between the perturbed TE and 
TM modes. For the perturbed TE modes we must have 

A™ 4- Bi 01 = -i(fi /Po)[Al 0) + B™] = 0, (105) 

A™ 4-£< 0) = coB(kUw), (106) 

while for the perturbed TM modes 

A r 4- Bl 0) = -i(p /Po)[A 3 0) 4- B™] = (KJKo) cos (kl Q w), (107) 
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AT + Bi 0) = 0. (108) 

If we now add to equation (103) the equivalent equation in B a , and 
make use of the fact that A™ + B™ t a = 1, 2, 3 are prescribed for the 
perturbed TE modes, in equations (105) and (106), we get a new equation 
involving only A { 2 0> and B™ . This equation together with equation (10G) 
can be solved for Ai 0) and B^ to yield (47) and (50). Once Ai 0) is known, 
A[ 0) and A< 0) can be determined from equations (103) and (48) yielding 
(46). We get B\ 0) and B l 3 0> from equation (105). In the same fashion we 
determine A™ and B™, a = 1, 2, 3 for the perturbed TM modes. 
Equations (102) and (104) [and the two equivalent equations in B[ l) 
and B^ l) ] are two inhomogeneous equations whose determinant vanishes. 
Thus the left side of (102) is a multiple of the left side of (104), and the 
equations are compatible only if the right side of (102) is the same 
multiple of the right side of (104) . This can be shown to be the case, and 
so (102) and (104) provide just one relationship between A[ l) and AJ°. 
There is a corresponding relationship between B[" and B"*. 

By replacing A a) B a , p, q by C a , D a , -r, -s, respectively, in the 
equations so far obtained, the formulas for the region £ ^ -a are ob- 
tained. Here — r and — s are the remaining two roots of the quartic 

Dip, j8) = 0. 

Next, if the assumed expressions for e a in |£| < a given by equation 
(11) are substituted into equations (90) through (92) we get four sets 
of three homogeneous, linear equations in F a , G a , L a , and M a , 
respectively, which hold for both the perturbed TE and TM modes. 
These equations are obtained from equations (95) through (97) by 
replacing A a and p by F a and -if, G a and ig, L a and - il, and M a and 
im, respectively. 

The determinental equation for each of these four sets of homo- 
geneous equations can again be expanded in powers of 77, and the coef- 
ficient of each power of rj separately equated to zero. The vanishing 
of the zeroth order coefficients yields equations (34), (41), (36), and 
(42) relating U,go,h, and ?» to O . The vanishing of the first order 
coefficients yields equations (44) and (45) relating /, , r/i , k , and m, 
to ft. 

Each of the four sets of homogeneous equations can be expanded in 
powers of 77, just as for the equations describing the region £ ^ a. To 
proceed further, we must again differentiate between the perturbed 
TE and TM modes. For the perturbed TE modes, equations (52) 
through (54) must be satisfied, while for the perturbed TM modes, 
equations (74) through (77) must be satisfied. These values satisfy the 
lowest order equations identically. 
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For both the perturbed TE and TM modes, the first order equations 
can be written as 

{K x - ffiFi" - /SjoFi" = -Sufi?, (109) 

-AWT + (K, - fl)F l 3 0) = -S 33 F?\ (110) 

(K 2 -11- fi)F? 

= -S ia F™ - S 23 F l 3 0> - (S 22 - 2/o/, - 2p P l )Fi°\ (111) 

For both the perturbed TE and TM modes, equation (111) vanishes and 
yields no information, while equations (109) and (110) have a nonzero 
determinant and so can be solved for F| n and Fi li yielding the solutions 
given in (55), (57), and (78). If we replace F a and / by G a and -g, we 
get equations which can be solved for G[ l) and G 3 l) yielding solutions 
given in (55), (57), and (78). The equations obtained when F a and / are 
replaced by L a and I, and M a and — m, respectively, have a different 
character. The two equations in L? and M?> corresponding to (111) do 
not vanish identically and can be solved for L' 1 ' and M { 2 \ The solutions 
are given in (61) and (80). The equations in L, (1) and L l 3 \ and M\ x) and 
M { 3 \ have a vanishing determinant. In the perturbed TE case, the 
equations are homogeneous and yield (59) and (60). In the perturbed 
TM case, the equations are nonhomogeneous but compatible, and yield 
the relations 

-hkLt" 4- (K 3 ti/KM l) 

= -*(& 3 - /o/3, - Z,j8 ) - }(&. - 2«aoAVft>ffs), (112) 

= -*(& 8 + / /3, + nhPo) + h(S 33 - 2l Q m i )(l K,/fi K 3 ). (113) 

We finally turn to the boundary conditions at £ = ±<r, of which there 
are eight, four at each boundary. They can be grouped as follows 

A 2 +B 2 = F 2 e ir ' + Q+*" + LJ" + il/ 2 <f '"", (114) 

C 2 4- D 2 = F 2 e~ if ' 4- Qj" + L 2 e~ ila + M 2 e i ""' , (115) 

-pA 2 - qB 2 = ifF 2 e if ° - igQ**" + ULJ 1 ' - imM 2 e ima ', (116) 

rC 2 4- sD 2 = ijF 2 e~ if ' - igG 2 e iaa + ilL 2 e~ iU - imM 2 e in ", (117) 

A 3 + B 3 = F 3 e ila + Q+'" + L 3 e { " + M 3 e im \ (118) 

C 3 4- D 3 = F 3 e~ i!tt 4- G 3 e<" + Lae'"' + Mae'", (119) 
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- V A, - qB 3 + 0(A 1 + BO = i(jF 3 + fiFJe"' + i(-gG 3 + fl?0<f *'* 
+ 2"(^3 + PLy + i(-mM 3 + /3il/ 1 )e~' mff , (120) 
rC 3 + sA, + z/S(C, + D,) - t'(/jFi + /3F0e _i/ " + ?(-<7G 3 + 0Gy° 

4- z(JL 3 + PLJe-"' + i(-mM s + /3il/,)e'" mff . (121) 
These equations split naturally into two groups, one group involving 
only the subscript 2 and the other group involving only the subscripts 
1 and 3. These equations can be expanded in powers of 77. The zeroth 
order equations are satisfied as long as (35) holds in the perturbed 
TE case and (37) holds in the perturbed TM case. 

For the perturbed TE modes, the first order expansion of equations 
(114) through (117) yields four nonhomogeneous equations in Ai 1 ' + 
B< n , C.<" + D< n , Fi li and G? } . The inhomogeneous terms on the'right 
side of these equations contain the parameter ft. The determinant of the 
equations vanishes, and then the condition that they be compatible 
provides an equation from which ft, given in (43), is determined. Once ft 
is determined, these equations yield (56) and (66). We can now choose 
the arbitrary parameter a y — indicated in (94) — so that F { 2 1) = 0. Then 
from (56) and (66) GJ" = A™ 4- 5<'* = 0. In addition, since ft is real, 
it can now be shown that r, = p*, s, = q*, C ( a 0) = A a n) *, D a 0) = B a 0) *, 
and C a X) + D a l) = [A l a u + B^]*, a = 1, 2, 3, which justifies equation 
(39). Finally, the first order expansion of equations (118) through (121) 
can be combined with equations (59) and (60), equation (102), and the 
corresponding three equations in B\ l) and B\ l) , C{ 1} and C 3 A) , and Z){" 
and D{" to form a set of equations from which A™ + B™ = [C™ 
+ Da*]*} and M" 1 , a = 1, 3, can be determined. These are listed in 
(58), (65), and (67). 

For the perturbed TM modes the procedure is virtually the same, 
except that it is now the first order expansion of equations (118) 
through (121) which has a vanishing determinant. The condition 
that these be compatible then yields the expressions (68) and (69) 
for fix . This set of equations also yields the result that 

L 3 l) + M<" = 0. (122) 

We can now pick the arbitrary parameter a t so that L^ = 0, which 
combined with (122) yields (85). Equations (112) and (113) now yield 
(86) and (87). The first order term of equation (118) then yields (88), and 
this result, combined with the equation obtained by adding equation 
(102) to the corresponding equation in B yields (89). Finally, equations 
(114) through (117) yield expressions (79) through (84) for F?\ G?\ 
L'", M< 1} , A?> + B™. 
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